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A hybrid algorithm is proposed for pure SU(N) lattice gauge theory based on Genetic Algorithms (GA)s and 
the Metropolis method. We apply the hybrid GA to pure SU(2) gauge theory on a 2-dimensional lattice and 
find the action per plaquette and Wilson loops being consistent with those given by the Metropolis and Heatbath 
methods. The thermalization of this newly proposed Hybrid GA is quite faster than in the Metropolis algorithm. 



1. Introduction 

The study of Genetic Algorithms based upon 
the theory of evolution originates with John Hol- 
land |lj in the mid-1970s. Applications in various 
fields have been proposed. In generic GAs, poten- 
tial solutions are represented as symbolic strings 
and operators are defined on the analogy of se- 
lection and genetic mechanisms in Nature. A 
crossover operator of exchanging genes between 
individuals leads to the special feature of GAs 
in which discrete points in searching space are 
treated at once, while normal Metropolis (MP) or 
heat bath (HB) algorithms treat only one point. 
We proposed a hybrid GA and applied it to pure 
SU(2) gauge theory in 2 dimensions. Physical 
values obtained by our hybrid GAs, such as the 
action per plaquette and Wilson loops at each 
(3 are consistent with those given by HBs. Our 
results show that the distinctive features of this 
new GAs leads to be the fast thermalization. The 
short calculation time are accomplished by the 
special encoding configurations on a lattice. 

2. Genetic Algorithms for SU(N) Lattice 
gauge theory 

In the language of GAs, a whole lattice corre- 
sponds to an individual, and the information on 
a field configuration is encoded in bit patterns 
which are treated as chromosome. The evalua- 
tion (fitness) function is the usual lattice action 
of SU(N) lattice gauge theory, 
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where U p is the element of SU(N) defined on a 
plaquette p and (3 is a coupling constant multi- 
plying S in the Boltzmann factor. 

We would like to call our newly proposed algo- 
rithm "Hybrid Genetic Algorithm" (HGA), be- 
cause it mixes a population update step of GAs 
with a normal Metropolis step applied to all mem- 
bers of the population. It combines heuristic 
methods for local search with GAs for global 
search with their respective advantages. In our 
algorithm, the GA starting from the t th genera- 
tion is followed by the Metropolis step, leading to 
the next (t + l) th generation, namely 
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where Pt denotes the population at t th genera- 
tion and Pgene pool represents a number of 2N pop 
genomes selected by Stochastic Universal Sam- 
pling which ensures unbiased selection ||. In 
this scheme, among all genomes of Pt laid on 
a pie chart with an area proportional to their 
fitness, 2N pop members are picked up simulta- 
neously into a gene pool by a single spin of a 
roulette wheel with 2 N pop equally-spaced point- 
ers. Two genomes taken from the gene pool are 
recombined by the crossover operator to give two 
offsprings (children). The children are subject to 
a mutation step which consists of bitwise random 
flips described by a (small) mutation rate. To 
establish thermal equilibrium, we set up an ac- 
cept/reject criterion in the next stage to fulfill 
the detailed balance condition: one does not re- 
place automatically old individual of the previous 
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Figure 1 . Thermalization of the action per plaquette on 32 • 32 lattice 



generation by new one but the better parent is 
replaced by the better child when the following 
condition is satisfied: 

min{l,exp~^ s ' !hUd ~ s ' aren ^} > £, (3) 

where £ is a uniform random number The 
population Pt+i is then the result of one normal 
Metropolis step. 

3. Thermalization 

In order to compare the thermalization rate 
with MP, a 2-point crossover scheme is adopted 
with a crossrate 0.65 and a mutation rate 0.008. 
The population sizes are 16 for j3 = 0.5, 64 for 
(3 = 2.0 and 128 for (3 = 8.0. Fig. [I] shows the 
results of these runs for three (3 values. In each 
figure, the horizontal dotted line shows the av- 
eraged action per plaquette for this (3 over the 
last 1000 sweeps of 30, 000 sweeps by the Heat- 
bath method. Square dots show the time his- 
tory corresponding to our HGA, the lines show 
the MC history of pure MP. The convergence 
to equilibrium (the target values given by HB) 
for the runs of our HGA is faster in CPU time 
than MP. In particular, for large (3 it is quite fast 
not being stuck due to low acceptance. The ap- 
proach to equilibrium in the run with = 0.5 
is, however, non-monotonous due to hard fluctu- 
ations. In the case of small /3s, most genomes 
in a population have very similar actions, which 
could have very similar genes or very different 
ones, an unexpected change of action is likely to 
be produced and accepted by the crossover opera- 
tion. To avoid strong fluctuations caused by bad 



crossing, a closer analysis of the correlation be- 
tween individuals is needed. Note that while MP 
treats just one lattice, HGA treats a population 
size number of lattices at once. 
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Figure 2. Comparison of Recombination 

schemes and population sizes at (3 = 2.0 

4. Recombination Schemes and Genetic 
Parameters 

We have tested two recombination schemes uni- 
form crossover (UC) and 2-point crossover (2C) 
with respect to their effectivity. UC in which ev- 
ery locus of a genome is randomly occupied by 
a gene from one of the parents with a probabil- 
ity p (or from the other with 1 — p), leads to an 
increased diversity of offsprings. In the lattice 
gauge theory it can lead to an increase of action, 
caused by a frequent mismatch of the links. In 
iVc-point crossover a new genome is created from 
partitions of parents genomes. It has less advan- 
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Figure 3. Physical values on 32 • 32 lattice 
tage for the diversity, but it reduces the mismatch 

problem of UC. To compare these schemes and 
to study the effect of population size for even- 
tual premature convergence, simulations are per- 
formed with f3 = 2.0, with two population sizes, 
32 and 16, with an equal crossrate of 0.65. Fig. |^ 
shows the comparison of schemes and population 
sizes for the case of f3 = 2.0. One sees that a 
large population size is useful for rapid decrease of 
the action (measured by the necessary number of 
generations), because it easier escapes from local 
minima thanks to the richer diversity in the popu- 
lation. On the point of elapsed time, however, the 
remarkable difference of convergence speed does 
not appear because of the large number of proce- 
dure. It suggests us the necessity of optimization 
of GA operators. Runs with UC did not converge 
fast, because of increasing link mismatch action 
that tends to outweigh the eventual decrease of 
bulk action. 



5. Conclusion 

We have designed a HGA and applied it to 2D 
SU (2) lattice gauge theory and calculated the ac- 
tion per plaquette and the Wilson loops. The 
results are shown in Fig. ||. For small (3 < 2.0 
the action per plaquette obtained by HGA have 
been averaged over the last 32 generations out of 
64 generations keeping a population size of 32, 
while for large /3s ( > 2.5) over 16 generations 
out of 32 generations with a population size of 
128. The average of Wilson loops obtained with 
MP were taken over 1000 iterations of the best 
genome among the population in the 50 th gener- 
ation of HGA for Ps between 2.0 < (3 < 4.0 (of the 
best genome in the 30 th generation for (3 > 4.5). 
Over the whole (3 region the result of the phys- 
ical averages are consistent with HB. They are 
obtained by HGA faster than by HB and nor- 
mal MP, without any further optimization of the 
GA part of our procedure. We have shown the 
possibility and effectiveness of GA schemes for 
SU(N) lattice gauge theory. A more detailed dis- 
cussion and analysis is needed to demonstrate the 
detailed balance properties. There is also enough 
room for further optimization of our algorithm. 
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